Bose-Einstein condensation of nonzero-center-of-mass-momentum Cooper pairs 



J. Batle, M. Casas 
Departament de Fisica, Universitat de les Illes Balears 
07071 Palma de Mallorca, Spain 

M. Fortes, M.A. Soli's 
Instituto de Fisica, Universidad Nacional Autonoma de Mexico 
Apdo. Postal 20-364, 01000 Mexico, DF, Mexico 

M. de Llano, A. A. Valladares 
Instituto de Investigaciones en Materiales, Universidad Nacional Autonoma de Mexico 
Apdo. Postal 70-360, 04510 Mexico, DF, Mexico 

O. Rojo 

PESTIC, Secretaria Academica & CINVESTAV - IPN, 04430 Mexico, DF, Mexico 

Received (January 15, 2001) 



Abstract 

Cooper pair (CP) binding with both zero and nonzero 
center-of-mass momenta (CMM) is studied with a set 
of renormalized equations assuming a short-ranged (at- 
tractive) pairwise interfermion interaction. Expanding 
the associated dispersion relation in 2D in powers of the 
CMM, in weak-to-moderate coupling a term linear in 
the CMM dominates the pair excitation energy, while 
the quadratic behavior usually assumed in Bose-Einstein 
(BE)-condensation studies prevails for any coupling only 
in the limit of zero Fermi velocity when the Fermi sea 
disappears, i.e., in vacuum. In 3D this same behavior is 
observed numerically. The linear term, moreover, exhibits 
CP breakup beyond a threshold CMM value which van- 
ishes with coupling. This makes all the excited (nonzero- 
CMM) BE levels with preformed CPs collapse into a single 
ground level so that a BCS condensate (where only zero 
CMM CPs are usually allowed) appears in zero coupling to 
be a special case in either 2D or 3D of the BE condensate 
of linear-dispersion-relation CPs. 
PACS: 74.20.Fg; 64.90+b; 05.30.Fk; 05.30.Jp 
Keywords: Cooper pairs; Bose-Einstein condensation; cu- 
prate superconductivity 

1. Introduction 

We consider an s-wave short-range, attractive 
(rank-one) separable interfermionic potential in d- 
dimensional momentum space Vpq = —{vQ/L'^)gpgq 
where uq ^ is the interaction strength, L the size 
of the system, and the ^p's are dimensionless form 
factors of the type = (1 +p'^ /pq)~^^^ in which po is 
the inverse range of the potential. Thus, e.g., po ~> oo 
implies g^ ~ 1 which corresponds to a contact or delta 



potential —v^biy) in configuration space. In either 2D 
or 3D such a potential well has an infinite number of 
bound states. As a result a many-fermion system with 
this interfermion interaction will collapse in the ther- 
modynamic limit to infinite binding per particle and 
infinite density. However, the potential can be "reg- 
ularized," i.e., constructed with infinitesimally 
small so that it supports a single bound state. 

The CP equation |^ for two interacting electrons 
of mass m above the Fermi surface, with momenta 
wavevectors ki and k2 and finite, nonzero CMM 
wavevector K = ki -I- k2, and relative momentum 
wavevector k = i(ki — k2), gives the total pair energy 
Ek = 2Ef — Ak in terms of vq, with Ep = h?kp/2m 
the Fermi energy. Here ^ is the CP binding en- 
ergy; it should not be confused with the BCS energy 
gap A. One can eliminate the variable vq in favor, in 
2D, of the vacuum bound-state energy i?2 5^ of the 
potential by combining Q the CP equation with the 
respective Lippmann-Schwinger one for the same in- 
terfermion interaction acting not in the Fermi sea but 
in vacuum. Then Ax can be extracted as a function 
of i?2 from the resulting renormalized CP equation 
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fPk^/m + Ak~ 2Ef + h?K^/4m 

2. Cooper pair dispersion relation 

After some algebra one finds the remarkable iden- 
tity, but only in 2D, that Aq = B2, i.e., for an attrac- 



tive delta interaction (regularized or not) the vacuum 
and zero-CMM CP binding energies coincide for all 
coupling. Using Ep/kp = hvp/2 one can expand Ak 
in powers of K for any coupling B2 and get 



(Ao - Ak) = -hvpK 
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h^K"^ /2{2m) will prevail for any nonzero coupling only 
when Ep/kp = hvp/2 0, i.e., in the vacuum limit 
when there is no Fermi sea. 

Figure 1 shows exact numerical results (full curves) 
of (1) in 2D for different B2/EP of the CP excitation 
energy ek/Aq as function of CMM K/kp. Note that 
the CPs break up at Ek/Aq = 1 where — 0, this 
being marked by large dots in the figure. In addition 
to the exact results we also exhibit the linear approxi- 
mation 2HvpK/tt (dot-dashed lines) for small B2/EF, 
as well as the quadratic approximation h?K^ /2{2m) 
(dashed parabolas) for large B2/EP. 

In 3D one obtains |^ similar results except that the 
dimensionless s-wave scattering length kpa'm. vacuum 
plays the role of a coupling parameter instead of the 
dimensionless binding energy B2/EP in the 2D case. 
Here, the limit Aq implies a ^ 0~ or 1/kpa 
—00 and corresponds to weak coupling, while the limit 
Aq — > 00 implies a — s- 0+ or l/kpa — > +00 and is 
strong coupling. In fact, for a = — |a| — > 0~ one finds 
Aq ^ (SEp/e"^) exp(— Tr/fcflol), a result first obtained 
by Van Hove On the other hand, a ^ 0+ yields 
Aq — > /ma^. Repeating the expansion carried out 
in 2D but without explicitly determining the coeffi- 
cient of the quadratic term gives 



ek = (Ao - Ak) 
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FIG. 1. Exact numerical results (full curves) for CP 
dispersion relation ek = Ao — Ak (in units of iC" = CP 
binding energy Ao) obtained from (|^) when gu = 1 for 
different coupling values B2/Ef- CPs break up when Ak 
turns negative, as indicated by large dots. The dot-dashed 
line is the linear approximation (virtually coincident with 
the exact curve for all B2 < O.lEp) while the quadratic 
approximation is shown dashed (see text for details). 

where a nonnegative CP excitation energy ek has 
been defined. It is this excitation energy that en- 
ters in the BE distribution function in determining the 
critical temperature in a picture of superconductivity 
as a BE condensation (BEC) of CPs. The leading 
term in (^ is linear in CMM, followed by a quadratic 
term. The latter is precisely the kinetic energy of 
what was originally the ordinary CP (and now is what 
is sometimes called a "local pair") — namely the fa- 
miliar nonrelativistic energy of the composite pair of 
mass 2m in vacuum. This dispersion relation has 
been the starting point for virtually all BEC studies 
of superconductivity (see, e.g., among others) . 

However, it is clear from (|^) that the quadratic term 



i.e., the same result cited in 1964 in Ref. for the 
BCS model interaction. The linear terms in both 
Eqs. (|) and (|) are identical [H for the BCS 
model interaction in weak coupling. In this case 
gi, = e{h^k^/2m - max[0, {Ep - Hujd)])0{Ep + huJD- 
ft^fc^/2m) where 9{x) is the Heaviside step function 
and ujd the Debye frequency. It becomes ijfc = 1 as 
hiUD — + 00 

3. Boson number 

Using a statistical model |l3j guaranteeing both 
thermal and chemical equilibrium in an ideal boson- 
fermion mixture, the number of bosons Nb{T) formed 
within the A'^-fermion system, valid at and below the 
BEC transition temperature Tc, is 

Nb{T)^^[N-No{T)] 

^^[1- (T/Tp) ln(l + e-/3{Ao(T)/2-M(T)})] 

where No{T) is the number of unpaired fermions, 
Ao(T) the appropriate finite-T generalization ||T^ of 
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the CP K = binding energy, /3 = l/fc^T, and the 
ideal Fermi gas chemical potential /i(r) in 2D is given 
exactly by 



^(T) = f3~^ ln(e 



1) Ep. 



(5) 



Figure 2 illustrates the zero CMM CP binding energy 
Ao(r) for three values of B2/fi{T). 
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FIG. 2. Temperature dependence of 2D zero-CMM 
CP binding energy Ao(r) vs. T/T^ for several couplings 
Bi/niT), where kgT^ = ^^. 

At T = (I) becomes 

Nb{0) = NAo{0)/4:Ef = NB2/4:Ef {B2 < 2Ef) 

= N/2 {B2 > 2Ep). (6) 

This should be compared with the BCS theory esti- 
mate (Ref. g, p. 128) 

NBiO)--iA/EF)^^^N{B2/EF), (7) 

where here A is the BCS T = energy gap, and the 
exact 2D result [0 A = \/2EfB2 was used in the 
last step. Since Nb < N/2, the estimate implies a 
breakdown for B2 > Ef/2 in the BCS case. 

4. Critical temperature 

Neglecting the background unpaired fermions and 
modeling the entire system as a pure boson gas of un- 
breakable CPs but with temperature-dependent bo- 
son number density nB{T) = Nb{T)/L'^, the explicit 
BEC Tc-formula for linear dispersion bosons in 2D 
[ p^ becomes an implicit one by allowing hb to be T- 
dependent, namely 



Ay/ShvF I 77^ 

ir = ^TTK- — yriByic, 
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This differs from the familiar BEC 3D formula Tc ~ 
S.Slh^n^J'^ /niBkB for quadratic-dispersion bosons. 
Both equations are special cases of the more general 

expression [|l5| of the form Tc oc Ug'^ for any space 
dimensionality d > and any boson dispersion rela- 
tion sk oc with s > 0. Solving (^) with (Q) and 
(1) for = self-consistently gives Tc/Tf vs. B2/EF 
as displayed in Figure 3 and compared with empirical 
values for cuprates that range pi] from 0.01 — 0.1. 
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FIG. 3. Critical BEC temperatures (full curves), for 
the pure unbreakable-boson gas, in units either of Tp or 
= jj,{T)/kB, compared with the BCS result (slanted 
dashed curve), vs. dimensionless coupling B2liJ.iT). Em- 
pirical cuprate data are taken from Ref. [n6|. 

Also shown in the figure are the BCS theory Tc's 
(see also Ref. [|l^) obtained by solving the single im- 
plicit equation 



dx Tp 

— tanh X 

X 2Tc 



In 
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where 7 is the Euler constant. Note that fc^Tc 
{e^ /'it)^/2B2Ep as coupling goes to zero, and also that 
2A/kBTc 27r/e^ ~ 3.53. 



5. BCS and BE condensates 

Finally, Figure 4 depicts in either 2D or 3D both 
condensates, the BCS one with its single K — pair- 
correlation state and the BE condensate |l^ with both 
(ground) K = and several (excited) K > CP states 
that form a "band" (shown in the figure as a discrete 
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spectrum for clarity) extending up to the breakup 
state Kq defined by Axg — 0. For perfectly linear dis- 
persion CPs, i.e., in 2D ek = — = ihvpK/n, 
the breakup CMM wavenumber is then just Kq = 
TTAo/2hvF. As this vanishes with coupling all the ex- 
cited boson levels collapse downwards and merge with 
the ground K = 
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FIG. 4. BCS pair condensate of Cooper correlations 
(CCs) and BE condensate of CPs, both below Tc, com- 
pared as explained in text, along with their respective 
normal states at T > Tc. Horizontal ellipsis indicate a 
fractional particle occupation which is macroscopic, or sig- 
nificant compared with unity. 

level, i.e., the "band" shrinks to the single ground 
level. Thus, for zero coupling the BCS condensate 
appears to be a special case of the BE condensate 
provided that the BCS CCs are essentially CPs, as is 
widely believed. 

6. Discussion 

Besides including the background unpaired fer- 
mions in the real mixture problem with our simple 



initial s-wave interfermion interaction, further refine- 
ments pending are: i) realistic Fermi surfaces; ii) Van 
Hove singularities or other means of accounting for 
periodic-crystalline effects; as well as the following in- 
teractions iii) the all-important c?-wave; iv) residual 
interbosonic ones; and v) the crucial CP-fermion inter- 
action vertex. It is precisely the latter ingredient that 
enabled T.D. Lee and coworkers and Tolmachev 
more generally, to link BCS and BEC through a re- 
lation whereby the BE condensate fraction is propor- 
tional to the (BCS-like) fermionic gap A(T) squared. 
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